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Abstract —Energy harvesting (EH), which explores renewable 
energy as a supplementary power source, is a promising 5G 
technology to support the huge energy demand of 
heterogeneous cellular networks (HCNs). However, the random 
arrival of renewable energy brings new challenges to network 
management. Through adjusting the distribution of traffic load 
in spatial domain, traffic shaping helps to balance the cell-level 
power demand and supply, enhancing the utilization of 
renewable energy. In this paper, we investigate the power 
saving performance of traffic shaping in an analytical way. 
Specifically, an energy-optimal traffic shaping scheme (EOTS) 
is devised for HCNs with EH, whereby the on-off state of the 
off-grid small cell and the amount of offloading traffic are 
adjusted dynamically, based on the statistic information of 
energy arrival and traffic load. Numerical results are given, 
which show that EOTS scheme can significantly reduce the 
on-grid power consumption for the daily traffic and solar 
energy profiles, compared with the greedy method where users 
are always offloaded to the off-grid small cell with priority. 

L Introduction 

To accommodate the exponential increasing traffic demand, 
small cells (SCs) are expected to be densely deployed overlaid 
with the conventional macro base stations (MBSs), introducing 
the heterogeneous cellular network (HCN) architecture. With 
high network capacity, HCNs will play an important role for 
5G evolution HI. However, the huge energy demand and costly 
deployment of SCs will cause heavy burdens to the network 
operators. To effectively reduce network energy consumption, 
dynamic network planning has been proposed and extensively 
investigated, whereby base stations can be dynamically turned 
off with traffic variations El-i). 

In addition to energy saving methods, energy harvesting 
(EH) technology is also introduced into HCNs, whereby new 
types of SCs equipped with EH devices (like solar panels or 
wind turbines) can exploit renewable energy for power 
source 0. In addition, together with mmWave wireless 
backhaul, EH can remove the wired connection of SCs, 
enabling flexible and cost-efficient network deployment. 
Whereas, the random arrival of renewable energy poses 
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significant challenges, and the design of efficient network 
management schemes is a key issue for HCNs with EH 0 . 

Plenty of works have been done to improve the utilization 
of harvested energy 0-0. Energy cooperation among BSs 
adjusts the renewable energy supply of each BS to match the 
traffic demand, which can be realized through dedicated 
power line connections between BSs 0 or the wireless 
energy transfer technology ||2l. However, it is costly to 
deploy dedicated wired connections among BSs, and the 
efficiency of wireless energy transfer is still very limited. BS 
cooperation can also be leveraged for traffic shaping, which 
adjusts the traffic distribution in spatial domain to match the 
energy status 0 0. Eor the single-tier homogeneous 
network jointly powered by EH and power grid, traffic 
shaping is realized by adjusting the cell size, with the 
objective to minimize the on-grid power consumption 0. In 
addition, the traffic shaping scheme based on dynamic 
programming is proposed for an HCN consisting of a 
conventional MBS and a SC jointly powered by the power 
grid and EH, where users are dynamically offloaded from the 
MBS to the SC according to the traffic and energy status 0. 

In this work, the power saving performance of traffic 
shaping is investigated in an analytical way. We consider an 
HCN covered by a conventional on-grid MBS, where one 
off-grid SC is further deployed to enhance the network 
capacity through EH. Traffic shaping, which offloads traffic 
from the MBS to the SC opportunistically based on the 
renewable energy status, can improve the utilization of 
harvested energy while guaranteeing the quality of service 
(QoS) (i.e., rate outage probability). Under the semi-dynamic 
traffic and energy models, we derive the closed-form on-grid 
power saving gain of traffic shaping, where the additional 
cost caused by user handovers is also considered. Based on 
the derived power saving gain, an energy-optimal traffic 
shaping (EOTS) scheme is proposed, whereby the SC is 
dynamically switched on/off and the amount of traffic 
offloaded to the SC is also adjusted, with the variation of 
energy arrival rate. Eor the daily traffic and solar energy 
profiles, numerical results show that the EOTS can 
significantly reduce the on-grid power consumption, 
compared with the greedy method where users are always 



Fig. 1: Illustration of a HCN with various energy sources. 


offloaded to the off-grid small cell with priority. The main 
contribution of this work is that the proposed EOTS scheme 
offers insights for real network operations, e.g., whether the 
SCs should be activated and how much traffic should be 
offloaded to the SCs. 

The rest of the paper is organized as follows. Section HI] 
introduces the system model, followed by the analysis of the 
power demand and supply in Section |In| Then, the EOTS 
scheme is proposed in Section |IV] based on the derived on- 
grid power saving gain. Numerical results are presented in 
Section |V] Einally, we conclude the paper in Section |Vl| 

II. System Model 
A. HCNs with Hybrid Energy Supply 

Consider a HCN whose basic coverage is guaranteed by a 
conventional on-grid MBS, and an off-grid SC is further 
deployed to enhance the network capacity with harvested 
energy, as shown in Eig. [T] Denote the coverage radius of 
MBS and SC as Dm and Dg, respectively. We model the 
user distribution as a non-homogeneous Poisson Point 
Process (PPP) in the spatial domain, whose density is /5m(f) 
(outside of the SC range) and Ps{t) (inside of the SC) at time 
t. Users within the SC can be either served by the MBS or 
offloaded to the SC for power saving, whereas users located 
outside of the SC range are all associated with the MBS. 
Thus, all users can be classified into three types: (1) MMUs, 
located outside of the SC and served by the MBS; (2) MSUs, 
located in the SC but served by the MBS; and (3) SSUs, 
offloaded to the SC. According to the properties of PPP im, 
the distributions of SSUs and MSUs within the SC should 
also follow PPP with density (p{t)ps{t) and (1 — (p{t))ps{t) 
respectively, if users located within the SC are offloaded to 
the SC randomly with probability (p{t) (i.e., random 

offloading scheme). Note that the offloading probability (p{t) 
reflects the amount of traffic offloaded to the SC, which 
should be adjusted according to the energy arrival rate. 

Assume the MBS and SC use orthogonal bandwidth to avoid 
the cross-tier interference. The bandwidths for the MBS and 
SC are denoted as Wm and 144 respectively, which can be 
totally or partially utilized according to the traffic demand. 


Denote by Wgs < Wg the utilized bandwidth at the SC, which 
is equally allocated to all SSUs. As for the MBS, Wmm is 
utilized to serve MMUs, while extra bandwidth lUms is shared 
by MSUs. To guarantee the QoS, tUmm and w^s are adapted 
to traffic and energy variations, satisfying uimm + tt'ms < Wm. 

Semi-dynamic traffic and energy models are adopted, and 
the time line is divided into T periods. In each period, the user 
density and energy harvesting rate are assumed to be static, but 
may change over periods. We solve the problem in a specific 
period t, and the subscript t is omitted in the following. 


B. Wireless Communication Model 

Denote by Pthi and Pts the transmit power of the MBS and 
SC, respectively. If user u is served by the MBS, the received 
signal to interference plus noise ratio (SINR) is as follows: 


7m = 


P'Tm'^u tin 


’/in 


( 1 ) 


Wm {e^ + i){a'^wuY 

where Wu denote the bandwidth allocated to user u, dm.u is 
the transmission distance, am is the path loss exponent of 
MBSs, /im,ii reflects Rayleigh fading (i.e., /im,u follows 
exponential distribution with mean 1), denotes the density 
of addictive Gaussian noise, and Omic'^Wu) equals to the 
interference among MBSs. Eor MMUs and MSUs, the 
achievable rates are as follows 

Tmm = l0g2(l + 7m), for MMUs, 

rVmm “r i 

rnM, = — -Tlog 2 (l+ 7 m), for MSUs, 


Kms + 1 


where and ATmm and Aims denote the number of residual 
MMUs and MSUs (except user u) , respectively. 

If user u is offloaded to the SC, the received SINR is 


_ Pts'^u tig^u hg^u 

Wg {da + 1) {a^Wu) ’ 

where dg.u represents the transmission distance, as is the 
path loss exponent of SCs, /is,„ reflects the effect of 
Rayleigh fading, and Os{a^Wu) equals to the interference. 
Since SSUs share Wgg equally, the achievable data rate of 
SSU u is given by: 

on 

= IC "1 1 + 7ss), (4) 

Ass + f 

where Kgg denotes the number of residential SSUs. 


III. Analysis of Power Supply and Consumption 
A. BS Power Consumption 

When a BS is working in active mode, its power 
consumption is given by nni 

Pbs = Po + (5) 

where Pq is the static power consumed by devices like the 
air conditioner, ^ reflects the power amplifier efficiency, Pp 
denotes the transmit power level, W is the system bandwidth 
w is the bandwidth utilized. If the BS is turned off (i.e., sleep 
mode), the power consumed Pg is negligible compared with 
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Fig. 2; Renewable energy arrival and consumption process. 


Pq, and thus Pg is approximated as zero. In this work, we 
consider a constant power level (i.e., Pp is a constant), while 
the power consumption is controlled by adjusting the utilized 
bandwidth w. 

B. Energy Queue Analysis 

We apply discrete energy model to describe the EH process, 
with a unit of energy denoted as E. Renewable energy arrival 
follows Poisson process with rate Ae, which is saved in the 
battery for future use. The battery capacity is considered to be 
sufficiently large (i.e., no battery overflow). Then, the battery 
status of the SC is illustrated in Fig. |2l with users offloaded to 
the SC opportunistically when the battery is not empty. Note 
that the renewable energy supply and consumption process 
can be modeled as an M/D/1 queue, with arrival rate Ae and 
service rate /te- The physical meaning of p,E is the Energy 
Consumption Rate of the SC, which depends on the utilized 
bandwidth Wss according to Eq. ©; 

Bb = ^ ^Pos + -^/3 sPts^ , (6) 

where the subscript “s” denotes “SC”. 

Next, we analyze the renewable energy queue. If ^ > 1, 
the energy queue is not stable and the queue length will go 
to infinity. In this case, the harvested energy is sufficient, and 
the SC can be always active for traffic offloading. Otherwise, 
the SC may be shut down due to empty battery, and 
meanwhile, the SSUs will go back to the MBS for service. 
The corresponding probability can be obtained according to 
queueing theory 02: 

Poff = l-—• (7) 

PE 

Notice that handover procedures are conducted at the 
moments of shutdown and reactivation, which brings 
additional signaling overhead and power consumption. The 

SC is shut down when the energy queue length transits from 
_ Ae , __ 

1 to 0, with frequency e 1121 . Considering the duality 

between operations of shutdown and reactivation, the 
additional handover power consumption Pho can be derived; 

Pho = 2C'hoPie“^/rE, (8) 

where Cho is in J denoting the energy consumption of one 
handover process, pi denotes the probability that the energy 


queue length is 1. For the M/D/1 queue, pi is given by pi = 

(l-^)(e™ -l) m. 

IV. Power Saving Gain of Traffic Offfoading 
Motivated by the power consumption model Eq. ©, the 
bandwidth required to satisfy the outage probability constraints 
should be obtained for power consumption analysis. In this 
section, the outage probabilities of different users are derived 
first, based on which we obtain the on-grid power consumption 
gain through traffic offloading. Furthermore, the EOTS scheme 
is proposed, which dynamically turns on/off the SC and adjusts 
the amount of traffic offloaded to the SC. 


A. Outage Probability Analysis 

Denote Gss and Gms the outage probabilities of SSUs and 
MSUs, respectively. We apply Lemmas 1 and 2 to simplify 
the outage probability constraints Gss < V and Gms < V- 
Lemma 1. If ^ 0 (high signal to noise ratio (SNR)) 
and —>• 0 (sufficient bandwidth), the service outage 

constraint of SSUs Gss = ]P{fss < F/th} < p is equivalent to 

t^ss^ss ^ ( 9 ) 


where Pth is the required data rate, Wss = ^ i^ the 

expected bandwidth allocated to each SSU, and Tss denotes 
the spectrum efficiency of cell edge users: 


-1 , Pts as + 2 r] \ 

Tss + 1 ) 2 W's Ps“V' 


( 10 ) 


Proof : Please refer to Appendix I 

Lemma 2. If ^ 0 (high SNR) and 0 

(sufficient bandwidth), the closed-form outage probability of 
the MSUs can be obtained by approximating all MSUs 
located at the SC. The outage probability constraint 
Gms = P{rms < Pth} < p is equivalent to 


^ms'^ms 


> Rth, 


( 11 ) 


where w^s 


__ 

l + {l — ip)ps'irDs‘‘ ’ 


= log2 1 + 


VPtu 


+ 1)D„ 


( 12 ) 


and Dms denotes the distance between the MBS and SC. I 
0 and 0 hold when the data rate 

requirement is relatively low compared with the spectrum 
resource, which is reasonable since BSs usually serve many 
users at the same time in practical systems. Eqs. dH) and (dH) 
suggest that the average data rate of the non-cell-edge users 
should be no smaller than Pth. As the proof of Lemma 2 is 
similar to that of Lemma 1, the details are omitted. We 
evaluate the accuracy of Lemmas 1 and 2 in Section |V] In 
addition, the outage probability constraint of the MMUs can 
also be simplified in the same way as Lemma 1, by 
approximating the distribution of MMUs as a PPP with 
density p'^ — — D^)/dJ^. The constraint 

Gmm = IPifinm < 7?th} < Is equivalent to 




> 7?th, 


(13) 












TABLE I: Simulation Parameters lITOl 


where Wmm = 




= log2 1 + 


and Tmm 
-Ppm 


are given by 
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B. On-Grid Power Consumption 

Based on the simplified outage probability, we derive the 
on-grid power consumptions when the SC is active or turned 
off, respectively. According to Eq. (|6]l, the bandwidth utilized 
at the SC Wss is constrained by the energy consumption rate; 

Wss = (peE - Pos) ■ (15) 

Ps^s 

Based on Lemma 1, we obtain the maximum offloading ratio: 

ip = 


Rth 

PsTtDs 


(16) 


Besides, according to Eq. ([n]l, the minimum bandwidth 
needed at the MBS to serve the MSUs is given as follows: 


= — (l + (1 - (p)psTrDs'^) . 


(17) 


However, the SC may be shut down due to empty battery, and 
meanwhile extra bandwidth Wmso is utilized at the MBS to 
serve the SSUs. Similar to Eq. ( fTTI ). we obtain w^so- 

Rth 


^mso — 


(1 + ipPsTTDs^) . 


(18) 


Therefore, the average bandwidth utilized at the MBS is given 
by Wmm + Wmsa + Pofft^mso, and the average on-grid power 
consumption is 

Ppm 


PmBS ~ Pom + Pn 


Wn 


(Wmm + W msa 4- PofftUniso) 4" Pho; 

(19) 

where poff and Pho are given by Eq. O and Eq. (|8]l when Ae < 
Pe, i e., the energy queue is stable. Otherwise, the renewable 
energy is sufficient, poff = 0 and Pho = 0. By combining 
Eqs. (fT5ll-(fT9]), the average on-grid power consumption can be 
derived when the SC is applied for offloading, with respect to 
the energy consumption rate of the SC /te- 

When the SC is turned off completely, ip = 0, and the 
bandwidth utilized at the MBS to serve the MSUs is equivalent 
to Wmsa + Wmso- Thus, the on-grid power consumption is 
Ppm 


PmBS ~ Pom + Pn 


Wn 


(tUinm 4“ tUmsa 4“ tUniso) • (20) 


Therefore, the on-grid power saving gain through offloading 
traffic from the MBS to the SC can be obtained; 

APmbs = Pmbs “ Pmbs 

P’pm 




"/ 5 mtUniso(l P*ho ( 21 ) 


A A pRF p 

where the first term APjJJgg represents the radio frequency 
power saved at the MBS. Based on Eqs. ( I15II181 I. we have 


A pRF 
MBS 


J (eeMeP — (CebPjs -I- ■ 
CbeAbE — (CbePjs + 


VFxnTxns 


Me < Ab 
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Parameter 

Value 

Parameter 

Value 

Dm 

lOOOm 

Ds 

300m 

POm 

130W 

POs 

56W 

PTm 

20W 

PTa 

6.3W 

fm 

4.7 

Pa 

2.6 

OLm 

3.5 

Ois 

4 

Wm 

lOMHz 

Ws 

lOMHz 

Rth 

lOOkbps 

V 

0.05 

dm 

1000 

ds 

2000 

a 

-105dBm/MHz 

-f^ms 

600m 


which (^ee = Eq. (l22l l reflects the 

conversion rate of harvested energy into on-grid power, i.e., 
the linear relationship between APj^gg and Ae. Substituting 
APmbs Pho in Eq. (ISTT l with Eqs. (l22l l and (|8]l, the 
power saving gain APmbs can be obtained. We summarize 
the relationship between APmbs and ^e in Proposition 1. 

Proposition 1. Let k = CeePos + , then APmbs 

has following properties: 

1) APmbs increases linearly with pE for /te < Ae; 

2) APmbs increases with p,E if k> SAeCeo ; 

3) APmbs decreases with pE for pE > Ae if k < (1 — 

^)AEC'ho ; 

4) If (1 - 2)AEC'ho < K < SAeCeo, APmbs is a concave 
function of ^ for //e > Ae, and the optimal condition 
is 




. ,±E 

Me 



Ae 

Me 



(23) 


Proof: Please refer to Appendix iBl I 

As 0 < Wss < Ws and ip < 1, pE should satisfy 

■^ < Mb < -^ j-Pos -h min |l, {pairD'^ + l)| /3 sPts| ■ (24) 


Eurthermore, the optimal energy consumption rate pE can be 
obtained by combining Proposition 1 and Eq. (l24ll . Eor 
example, the SC should be active and work at the maximal 
energy consumption rate when Cho = 0 (i.e., no handover 
cost), according to the second property of Proposition 1. In 
addition, the SC should be activated only when the power 
saving gain is positive; otherwise, the activation will increase 
the on-grid power consumption. 

At this point, the energy-optimal traffic shaping (EOTS) 
scheme has been obtained, whereby the on-off state of the 
SC and the amount of traffic offloaded should be 
dynamically adjusted with the energy arrival rate, based on 
Proposition 1 and Eq. (l24ll . 


V. Simulation Results 

In this section, the accuracy of the derived outage 
probability in Lemmas 1 and 2 is validated, and the power 
saving performance of EOTS is also evaluated. Simulation 
parameters are listed in Table U Solar powered SC is 
considered. Eig. 0 shows the daily traffic ifTOl and solar 
energy profiles ifTSll . 
























Fig. 3: Daily traffic and energy profiles. 
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Fig. 4: Outage probability. 


A. Outage Probability Evaluation 

We evaluate the analytical results in Lemmas 1 and 2 via 
Monte Carlo simulations, with user density set as 
Pni=20/km^ and ps=60/km^. The number of users and their 
locations are generated randomly, and the results of 5000 
simulation samples are averaged. The analytical results 
obtained by Eqs. (© and (fTSl) are compared with the 
simulation ones in Fig. EJa), which are shown to increase 
with the data rate requirement. Fig. |4fb) validates Lemma 2 
with different distances between the MBS and the SC, when 



Fig. 5: Maximal power saving gain of traffic offloading. 



Fig. 6; Average power saving gain of EOTS. 

the bandwidth allocated to MSUs is set as Wms = 1 MHz. It 
can be seen that the analytical and simulation results match 
well, with acceptable deviations. 

B. Power Saving Gain through Traffic Shaping 

Fig.E] shows the power saving gain achieved by activating 
the SC for traffic offloading, with the optimal energy 
consumption rate. Eirstly, the power saving gain decreases 
with the handover cost Cho, which is straightforward. 
Secondly, the power saving gain firstly decreases and then 
increases with the energy arrival rate. The reason is that the 
on-grid power consumption is dominated by the handover 
cost when the energy arrival rate Ae is small, due to the 
insufficient green energy supply. Besides, the handover 
frequency increases with the energy arrival rate, given that 
the energy queue is empty during most of the time with low 
energy arrival rate. Accordingly, the total power saving gain 
decreases when Ae is small. Then, as Ae further increases, 
the harvested energy balances out the handover cost, which 
improves the power saving gain. 

Fig. |6] shows the power saving performance of EOTS, 
which is compared with the greedy scheme for the traffic 
and energy arrival profiles given by Eig. [3 Eor the greedy 
scheme, the SC is always active and work at the maximal 
energy consumption rate. The average power saving gain 
with respect to the maximum energy arrival rate is shown. 

































where the maximal user density is set as = 100/km^. 

When the handover cost is ignored (Cho = 0), the two 
schemes have the same performance. However, EOTS is 
more advantageous when the handover cost increases. Notice 
that, the power saving gain of the greedy scheme may even 
become negative, i.e., deploying a off-grid SC may increase 
the total power consumption without effective traffic 
management. Whereas, EOTS always reduces the on-grid 
power consumption. 

VI. Conclusions and Euture Work 


Eurthermore, > 0 (i.e., sufficient bandwidth), 2 — 1 = 

In 2 ■ Hence, Lemma 1 is proved. 

Appendix B 

Proof of Proposition 1 

Denote Q = ^ where g G (0,1) for /te > Ae, and we take 
the derivation of APmbs with respect to g: 


(IAPmbs 

dg 


= —K — AsCho 




= -« + AECho/(i?). 


In this paper, we derived the closed-form expression of the 
on-grid power saving gain achieved by traffic offloading, 
whereby users are dynamically offloaded from the 
conventional on-grid MBS to the off-grid SC powered solely 
by EH. The derived power saving gain indicates the 
conversion rate of harvested energy into on-grid power, with 
spatial traffic shaping. Besides, an energy-optimal traffic 
shaping scheme, EOTS, is proposed, which offers a 
guideline for practical energy-aware network operations. In 
the future, the co-existence of multiple SCs with diverse 
power sources will be studied. 

Appendix A 
Proof of Lemma 1 

Due to channel fading or bandwidth limitation, the 
achievable rate of a user may be smaller than its required 
threshold i?th, and correspondingly, this rate outage 
probability is defined as Gss = ]P{rss < ^th}- Next, we 
derive the closed-form expression of Gss, based on the 
Eqs. (O and (|4|i. As the distribution of SSUs follows PPP 
with density (pps, the probability distribution function of the 
distance between SSUs and the SC ds is fdsid) = Based 
on Eq. ©, we have 


■{ 


P i 7ss > 

rD, 


> - 1 


> (9. + lWWs /2(^-+i)5S _ 1 
- V 

(6>s -p lIo-^Ws A(/Ces + 1) 


2d 


-dd 


PTarf““= 
2t 


2d 


- 1 ) ) —dd (25a) 


^ ^ {es + l)a^Ws _ A 

Os -I- 2 Pts V / 


where Eq. (125 al l is due to the property of Rayleigh fading, 
and (I25bl i comes from 0 (i.e., high signal to noise 

ratio). In addition. Kgs should follow the Poisson distribution 
of parameter nippsDs'^, according to the Slivnyak-Mecke 
theorem CD. Therefore, the outage probability of SSU Gss 
satisfies 


Gss = l- ^ P (^7bb > 2^^+^^^ - Pk,AI<) 
= 1- £ F (.ss > 2^-+^)^ - l) 

ir—n A / 


(26) 


2P>A(6»s + 1)(72Ws 


Pts (os-I-2) 


2 exp ( 7rDg ps ( 2 “ss — 1 ) ) — 1 


Notice that 

dp ,3 [ 2 


(28a) 


we have 1 — e ^ < f{g) < 3/2 for 0 < p < 1. In addition, 
APmbs is concave with respect to g, as = f'id) < 

0. Specifically, < 0 for k > SAeGeo, and > 

0 for K < (1 — i)AEGho. Otherwise, there exists g satisfying 
^A|kBS.|~ = 0, and the corresponding energy consumption 
rate ^ can maximize APmbs- Hence, Proposition 1 is 
proved. 
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